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FOREWORD 


This  final  research  report  was  prepared  jointly  by  the  School  of 
Electrical  Engineering  and  the  Engineering  Experiment  Station  of  the 
Georgia  Institute  of  Technology,  Atlanta,  Georgia,  under  Contract  No. 
DAAH01-73-C-0796 ,  This  contract  was  initiated  by  the  RF  Guidance  Tech¬ 
nology  Branch  of  the  Advanced  Sensors  Directorate,  Research,  Development, 
Engineering  and  Missile  Systems  Laboratory,  U.  S.  Army  Missile  Command, 
Redstone  Arsenal,  Alabama.  The  contract  was  administered  under  the 
direction  of  Carlton  H.  Cash  and  Thomas  W.  Morgan  of  the  RF  Guidance 
Technology  Branch. 

The  period  of  performance  covered  by  this  report  extends  from 
1  May  1973  to  31  December  1973. 

Report  authors  are:  Edward  B.  Joy  and  Gene  K.  Huddleston  of  the 
School  of  Electrical  Engineering  and  Harold  L.  Bassett,  Charles  W.  Gorton 
and  Steve  H.  Bomar,  Jr.,  of  the  Engineering  Experiment  Station.  Other 
personnel  participating  in  this  work  includes  Earle  A.  Welsh,  Research 
Engineer,  and  James  N.  Farmer,  Student  Assistant. 

The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  necessarily  representing 
the  official  policies,  either  expressed  or  implied,  of  the  U.  S.  Army 
Missile  Command  or  the  U.  S.  Government. 
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SECTION  I 


INTRODUCTION 

This  report  describes  the  results  of  a  research  program  to  select, 
analyze,  and  evaluate  various  radome  materials  and  configurations  in  order 
to  determine  a  preliminary  design  of  a  radome  which  can  be  usee  with  advanced 
RF  seekers  operating  over  the  8.0  to  18.0  GHz  frequency  band  and  physically 
constrained  for  operation  in  a  five-inch  (outside)  diameter  missile. 

The  basic  radome  configuration  considered  consists  of  a  f  .agent  ogive 

shape  surmounting  a  cylinder  and  an  n-order  monolithic  wall  structure. 

«• 

Attention  ^s  directed  toward  the  first-order  (half-wave)  wall  structure 
because  of  it'-  superior  bandwidth  compared  to  the  higher  order  structures; 
furthermore,  the  5-inch  outside  diameter  requirement  dictates  as  thin  a  wall 
as  possible  to  allow  room  inside  for  the  antenna  system.  Details  of  the 
radome  shape  and  position  of  the  antenna  inside  the  radome  are  given  in  the 
drawing  of  Figure  11. 

The  radome  design  process  described  herein  follows  the  natural  course 
of  determining  the  temperatures  and  mechanical  stresses  to  be  encountered, 
selecting  candidate  radome  materials,  and  assessing  the  effects  of  radome 
chape  and  radome  material  on  the  electrical  performance  of  the  seeker  system. 

■1  4 

Based  on  these  results,  a  radome  design  is  selected  which  yields  the  best 
electrical  performance  while  meeting  the  requirements  of  the  flight  environ¬ 
ment  . 

Section  II  presents  a  discussion  of  the  assumed  aerodynamic  heating 
and  drag  considerations  for  the  radome  under  study.  The  range  of  temperatures 


1 


and  effects  of  radome  fineness  ratio  on  missile  trajectory  and  range  are 
also  presented . 

Section  III  presents  electrical  and  mechanical  data  for  candidate 
radome  materials.  Criteria  for  selecting  a  single  material  from  among 
those  studied  are  also  discussed . 

Section  IV  presents  computed  electrical  performance  data  for  a  number 
of  radome  configurations . 

The  conclusions  drawn  from  the  data  presented  are  given  in  Section  V. 
A  radome  design  is  also  recommended. 


9 


liHi'.r.'i.lmrf 


SECTION  II 


AERODYNAMIC  l -EATING  AND  DRAG 

1.  Introduction 

One  of  the  first  steps  in  the  design  of  any  radome  is  the  determina¬ 
tion  of  the  effects  of  the  flight  environment  on  the  radome.  One  of  the 
most  important  effects  to  be  determined  is  the  range  of  temperatures  and 
stresses  to  be  encountered  by  the  radome  following  a  specified  trajectory. 
Another  important  effect  is  that  of  the  shape  of  the  radome  on  the  missile 
trajectory.  These  effects  of  temperature,  stress,  and  radome  shape  are 
important  because  they  influence  greatly  the  choice  of  radome  material 
and,  ultimately,  the  electrical  performance  of  the  seeker  antenna  system 
enclosed  by  the  radome. 

In  this  section,  the  aerodynamic  heating  and  drag  considerations  are 
examined.  Upper  bounds  on  radome  temperatures  are  calculated  independently 
of  the  radome  material,  radome  shape,  and  trajectory  by  using  a  steady- 
state  analysis  and  a  flat  plate  assumption.  These  results,  presented  in 
Table  I,  influence  the  initial  choice  of  candidate  radome  materials. 
Following  the  calculation  of  upper  bounds  on  the  temperature,  the  results 
of  a  transient  temperature  analysis  are  presented  for  three  candidate 
materials.  The  transient  analysis  does  account  for  the  flight  trajectory 
and  radome  material.  The  effects  of  the  radome  shape  and  wall  thickness 
on  the  mechanical  stresses  and  missile  trajectory  are  also  calculated  and 
presented . 
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The  initial  calculations  related  to  aerodynamic  heating  were  the  deter¬ 
mination  of  the  stagnation  temperature,  the  flat  plate  laminar  recovery 
temperature,  and  the  flat  plate  turbulent  recovery  temperature.  The  stagna¬ 
tion  temperature  is  the  upper  bound  on  the  temperature  that  any  location 
on  the  radome  surface  can  reach.  Physically,  it  is  the  temperature  attained 
by  the  air  flowing  around  the  radome  (at  the  outer  edge  of  the  boundary 
layer)  at  the  stagnation  point  on  the  radome.  For  a  radome  at  zero  angle 
of  attack,  the  stagnation  point  is  the  radome  tip.  Actually,  the  outer 
surface  of  the  radome  at  the  stagnation  point  will  reach  a  temperature 
that  is  les9  than  the  stagnation  temperature  because  of  thermal  radiation 
from  the  surface  and  heat  conduction  into  the  radome  material.  The  recovery 
temperatures  calculated  correspond  to  the  temperature  of  a  thin  plate 
placed  parallel  to  the  flow  under  adiabatic  conditions;  i.e.,  internal  con¬ 
duction  and  thermal  radiation  are  assumed  to  be  zero.  Physically,  the 
recovery  temperatures  represent  an  upper  bound  of  the  outer  surface  tempera¬ 
ture  in  the  aft  region  of  a  pointed  radome.  Whether  the  laminar  or  turbulent 
recovery  temperature  is  considered  depends  on  whether  the  boundary  layer 
is  laminar  or  turbulent  at  the  point  under  consideration.  The  results  of 
the  calculations  of  stagnation  and  recovery  temperatures  are  presented  in 
Table  I.  These  calculations  are  based  on  sea-level  conditions  for  a 
standard  day  under  steady-state  conditions. 

The  transient  temperature  response  of  the  hemispherical  radome  was 
determined  at  the  stagnation  point  by  means  of  finite  difference  equations 
and  the  aid  of  a  high-speed  digital  computer.  Inputs  to  the  computer-aided 


Flight  Velocity 
(meters/sec.) 

Stagnation 
(  F) 

Temperatures 

Laminar  Recovery 
(°F) 

Turbulent  Recovery 
- - 

0 

59 

59 

59 

200 

95 

90 

91 

400 

202 

180 

185 

600 

380 

330 

342 

800 

623 

537 

555 

1000 

925 

796 

826 

temperature  analysis  included  radome  material  properties,  nose  radius, 
radome  wall  thickness,  and  flight  velocity  as  a  function  of  time.  The 
radome  material  properties  needed  for  determining  temperature  profiles  are 
density,  specific  heat,  thermal  conductivity,  and  surface  emittance.  In 
addition,  the  coefficient  of  thermal  expansion,  modulus  of  elasticity,  and 
Poisson's  ratio  were  used  as  inputs  so  that  thermal  stresses  could  be  cal¬ 
culated.  Also,  a  heat  transfer  coefficient  for  the  outside  surface  (front 
side)  arid  a  boundary  condition  for  the  inside  surface  (back  side)  of  the 
radome  were  specified.  The  heat  transfer  coefficient  was  calculated  using 
an  equation  relating  nose  radius  and  flight  velocity  for  sea  level  flight. 
The  equation  used  is  an  empirical  relationship  which  was  developed  following 
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the  procedure  recommended  by  van  Driest  [1].  The  back  side  of  the  radome 
was  assumed  to  be  perfectly  insulated;  this  is  a  fairly  universal  assump¬ 
tion  which  uncouples  the  transient  radome  heating  from  the  thermal  response 
* 

of  the  contents  of  the  radome.  Although  a  variety  of  intermediate  results 
are  obtained,  those  of  primary  interest  are  front  side  temperature  (to  see 
if  material  temperature  limits  are  exceeded),  back  side  temperature  (to  see 
if  the  interior  of  the  radome  becomes  too  high  in  temperature) ,  and  thermal 
stress  (to  see  if  material  failure  will  result) .  The  calculations  for  the 
three  materials  of  interest  are  presented  in  Figures  1  and  2.  It  is  noted 
that  the  temperatures  and  thermal  stresses  shovm  in  these  figures  are 
well  within  the  operating  limits  of  the  three  materials  shown. 

3.  Radome  Drag  Coefficient 

As  the  fineness  ratio  L/D  of  the  radome  is  changed,  the  overall  drag 
of  the  hypothetical  missile  changes,  resulting  in  a  modified  trajectory. 

In  investigating  this  effect,  the  modified  trajectory  was  related  to  the 
original  trajectory.  The  assumption  was  made  that  the  overall  drag  co¬ 
efficient  was  constant  throughout  the  flight.  The  trajectory  calculations 
for  each  fineness  ratio  were  made  using  the  same  missile  weight,  fuel 
weight,  thrust,  and  burn  time.  It  is  noted  that  the  minimum  velocity 
required  for  maneuverability  is  reached  for  a  fineness  ratio  of  0.5  at 
about  0.67  on  the  normalized  time  axis  of  Figure  3.  The  results  of  these 
calculations  are  presented  in  Figure  3. 

In  order  to  evaluate  the  mechanical  loading  on  the  radome  with  an  angle 
of  attack,  calculations  were  made  using  a  ten  degree  angle  of  attack  at  the 
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maximum  velocity  in  the  trajectory.  Sea  level  conditions  were  assumed. 
These  calculations  are  based  on  a  Newtonian  pressure  distribution  using  the 
results  presented  by  Mayo  [2],  The  stresses  are  calculated  from  the  aero¬ 
dynamic  axial  load  and  bending  moment.  The  results  of  these  calculations 
are  presented  in  Table  II.  The  stresses  indicated  refer  to  the  maximum 
compressive  or  tensile  stress  in  the  radane  wall  at  the  base  of  the  radome. 
All  of  the  stresses  shown  in  Table  II  are  small  for  the  three  materials  of 
interest. 


TABLE  II 

MECHANICAL  LOADING  AT  MAXIMUM  VELOCITY 
,  ANGLE  OF  ATTACK  =  10° 


L/D 

Radome  Wall  Thickness 
( inches) 

Compressive  Stress 
(p.s.i.) 

Tensile  Stress 
(p.s.i.) 

% 

0.250 

224 

0 

L 

o 

0.375 

153 

0 

l 

0.250 

242 

<  1 

l 

0.375 

170 

<  1 

2 

0.250 

450 

348 

2 

0.375 

330 

242 

3 

0.250 

698 

642 

3 

0.375 

500 

462 
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SECTION  III 


RADOME  MATERIALS  SELECTION 

1.  Introduction 

The  thermal  and  mechanical  data  presented  in  Section  II  were  used  as 

t 

guidelines  in  selecting  candidate  radome  materials.  The  criteria  used  in 
selecting  the  materials  were:  (1)  strength  during  peak  temperatures,  (2) 
thermal  shock  resistance,  (3)  microwave  transmission  properties,  (4)  weight, 
(5)  fabricability,  (6)  rain  erosion  resistance,  and  (7)  cost.  Some  of  the 
materials  considered  are  listed  in  Table  III.  The  electrical  properties  of 
all  the  listed  materials  are  relatively  good  up  to  the  temperatures  indicated 
and  higher;  e.g.,  measurements  at  Georgia  Tech  on  polyimide-quartz  samples 
indicate  that  the  material  will  operate  to  at  le-'-'t  1000°F  with  only  small 
degradation  of  electrical  properties. 

For  background  information  and  to  put  some  of  these  materials  into 
perspective,  a  review  of  the  current  operating  radomes  is  in  order.  Glass 
reir  „d  plastic  radomes  have  been  widely  used  to  the  range  of  Mach  3. 

example,  the  B-l  aircraft  radome  is  constructed  of  polyimide  quartz, 
er  aircraft  have  used  radomes  made  of  glass  reinforced  polyesters, 
epoxies,  and  silicones.  Available  fibers  include  E-glass  (most  common), 

S -glass  (slightly  higher  temperature  capability  than  E-glass),  fused  quartz 
(.highest  temperature  capability  and  most  expensive),  and  PRD-49  (a  high 
temperature  organic  fiber  related  to  nylon  and  manufactured  by  DuPont). 

Other  resins  are  available  including  polyurethane  (low  temperature  capability 
used  on  commercial  aircraft),  silicones  (good  elevated  temperature  properties 
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silicon  nitride 


but  low  strength),  phenyl  silanes  (similar  fco  silicones,  used  on  the  F-14 
aircraft).  Teflon  (non-charring  ablator),  and  polybutadiene  (a  newly  developed 
high  temperature  resin,  manufactured  by  Firestone).  With  present  technology, 
none  of  the  multitude  of  fiber  reinforced  resin  systems  can  be  expected  to 
find  applications  above  Mach  3  to  3.5.  However,  this  class  of  materials  is 
experiencing  st  »dy  development  for  uses  up  to  this  range  of  Mach  numbers; 
polybutadiene  resin  and  PRD-49  are  the  newest  additions,  just  now  becoming 
available  in  commercial  quantities.  As  a  class,  the  glass  fiber  reinforced 
plastics  are  usually  preferred  over  ceramics  for  radome  applications  where 
mission  requirements  are  mild  enough  for  them  to  operate  successfully. 

They  do  not  suffer  from  the  brittleness  of  ceramics,  are  light  weight, 
relatively  easy  to  fabricate,  and  relatively  low  in  cost. 

Aluminum  oxide  has  been  used  for  missile  radomes  operating  up  to  about 
Mach  3.  One  example  is  the  Sparrow  III  air-to-air  missile.  This  material 
has  excellent  rain  erosion  resistci.ee,  but  above  Mach  3  it  cannot  survive 
the  thermal  shock  resulting  from  aerodynamic  heating. 

Cordierite  (PYUOCERAM  9606®  and  Raycerara  III)  has  been  used  for  radomes 
up  to  about  Mach  4.  It  is  used  on  the  Navy's  Standard  Missile  for  surface 
to  air,  anti-aircraft  operation.  Cordierite  has  good  rain  erosion  resistance, 
but  a  low  melting  point  and  moderate  thermal  shock  capability.  It  was  a 
leading  candidate  for  the  Army's  Sam-D  missile,  but  was  eventually  discarded 
in  favor  of  slip-cast  fused  silica  because  of  thermal  shock  considerations 
and  the  fact  that  it  is  manufactured  by  a  proprietary  process.  It  was  a 
candidate,  however,  for  this  program  and  was  carefully  considered  as  the 


radome  material. 


Slip-cast  fused  silica  presently  appears  to  be  the  only  candidate  radome 
material  for  use  above  Mach  6.  That  does  not  mean,  however,  that  it  is  ex¬ 
cluded  as  a  candidate  material  at  lower  speeds.  This  material  is  being 
used  on  the  SAM-D  program  and  has  been  a  candidate  for  other  missile  pro¬ 
grams.  This  material  has  excellent  thermal  shock  rtsistance,  but  has  poor 
rain  erosion  resistance.  Fabrication  technology  has  reached  the  point  of 
commercial  reliability,  and  the  cost  is  moderate. 

Boron  nitrides,  both  hot  pressed  and  isotropic  pyrolytic  boron  nitride, 
have  excellent  thermal  shock  resistance  but  poor  rain  erosion  resistance. 

They  hive  not  been  seriously  considered  for  radomes  because  they  cannot  be 
fabricated  in  large  shapes;  even  if  fabrication  technology  were  developed, 
the  cost  would  probably  be  prohibitive.  They  are  both  good  candidates  for 
small  windows  operating  to  Mach  8  or  higher. 

Beryllium  oxide  has  excellent  rain  erosion  resistance  and  good  thermal 
shock  resistance.  It  has  never  been  used  for  an  operational  radome,  pri¬ 
marily  because  of  the  toxic  nature  of  finely  divided  beryllium  oxide  powders. 
Costs  of  machining  with  the  required  safety  precautions  are  prohibitive. 

It  has  been  used  for  high  power  transmission  windows. 

Reaction  sintered  silicon  nitride  is  a  new  material  currently  under 
development  for  possible  radome  use.  Potentially,  it  offers  good  thermal 
shock  resistance  in  aerodynamic  heating,  much  better  rain  erosion  resis¬ 
tance  and  higher  strength  than  fused  silica.  It  probably  will  operate  into 
the  Mach  6  range.  Fabrication  techniques  and  cost  should  be  comparable  to 
fused  silica,  but  this  material  is  not  yet  ready  for  commercial  production. 
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From  the  preceding  discussion  of  candidate  radome  materials,  several 
observations  are  apparent: 

(1)  Various  materials,  including  glass  fiber  reinforced  plastics 
and  aluminum  oxide,  are  suitable  for  radomes  up  to  Mach  3. 

(2)  Cordierite  is  suitable  for  radomes  up  to  about  Mach  4  or  5- 

(3)  Slip-cast  fused  silica  is  presently  the  only  strong  candidate 
for  radomes  in  the  range  above  Mach  5.  This  material  has  re¬ 
latively  low  strength  and  poor  rain  erosion  resistance.  Reaction 
sintered  silicon  nitride  may  offer  relief  from  these  deficiencies, 
up  to  about  Mach  6,  with  further  development. 

(4)  Beryllium  oxide  and  the  boron  nitrides  are  not  suitable  for 
radomes  because  of  fabrication  difficulties  and  cost;  they  are 
useable  for  dielectric  windows  with  isotropic  pyrolytic  boron 
nitride  being  useable  up  to  very  high  Mach  numbers. 

2.  Materials  Data 

For  the  application  of  concern  here,  four  candidate  radome  materials 
were  considered:  polyimide-quartz,  glass  reinforced  polybenzimidazole, 
slip-cast  fused  silica  and  cordierite  (Pyroceram  9606^).  For  comparison 
purposes,  the  strength  properties  of  four  glass  fiber  reinforced  plastic 
materials  are  plotted  in  Figure  4.  The  strength  properties  of  the  PBI 
laminates  are  plotted  in  Figures  5  and  6.  The  ceramic  material  strength 
properties  are  plotted  in  Figure  7.  Alumina  is  included  here  for  com¬ 
parison  purposes.  The  radome  stresses  shown  in  Figure  2  of  Section  II  are 
very  low  compared  to  the  maximum  stress  that  the.  materials  will  withstand. 
Both  the  poJyimide/E-Glass  BPI-373  and  the  polyiraide  quartz  BPI-373  maintain 

i 

good  strength  properties  to  800°F.  The  PBI  laminate  possesses  slightly 
better  properties  as  indicated  in  Figure  5. 
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The  other  considerations  for  material  properties  were  the  electrical 
characteristics.  As  indicated  in  Figure  8,  the  dielectric  constant  of 
polyioide -quartz  BPI-373  increases  slightly  from  3.25  to  3.3  over  a 
temperat;  -e  range  from  ambient  to  600°F.  This  is  a  1.5  percent  increase 
and  will  net  cause  serious  degradations  in  radome-antenna  performance  due 
to  this  effect.  This  change  is  also  typical  of  PBI.  The  electrical  pro¬ 
perties  of  slip-cast  fused  silica  and  Pyroceram®  shown  in  Figures  9  and  10 
change  very  little  over  the  operating  temperature  range  shown  in  Figure  1. 
This  particular  property  is  a  feature  of  the  ceramic  materials,  and  it  is  an 
advantage  that  is  given  due  consideration  when  choosing  a  radome  material. 
Fortunately,  the  maximum  expected  radome  outer  surface  temperature  is  not 
very  high  so  that  the  organic  materials  considered  can  survive  the  environ¬ 
ment  . 

In  general,  the  expected  change  in  electrical  performance  of  the 
radome  due  to  the  temperature  dependence  of  the  electrical  properties  is 
not  significant  for  the  four  radome  materials  considered.  The  perused  data 
are  not  included  in  this  report  on  organic  materials  whose  temperature- 
electrical  performance  fell  below  the  two  organics  considered.  For  example, 
Epoxy/E-Glass  E-100  has  a  7  percent  variation  in  dielectric  constant  from 
ambient  to  600°F  and  its  loss  tangent  is  greater  than  0.03.  Similar  re¬ 
sults  have  been  found  with  Polyester /E-Glass,  Vibrin-135  and  Polyimide/ 
E-Glass  BPI-3  73  . 

Since  the  environment  is  such  that  an  organic  material  will  survive, 
it  is  advantageous  to  eliminate  the  ceramics  and  consider  only  the  organic 
materials  polyimide  quartz  and  PBl.‘  The  ceramics  are  much  more  fragile  than 
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Figure  10.  Dielectric  Constant  and  Loss  Tangent  of  Pyrocer am®  9606 
versus  Temperature  [3] 


the  organics;  and  if  it  is  found  that  grinding  of  the  ceramic  is  necessary 
to  obtain  a  desired  thickness,  the  organic  design  would  be  less  expensive. 
The  most  desirable  feature  of  Cordierite  in  its  resistance  to  rain  erosion. 
It  will  withstand  rainfall  conditions  much  better  than  the  other  materials. 
Slip-cast  fused  silica  will  absorb  moisture  and,  although  there  are  tech¬ 
niques  to  seal  the  radome,  this  adds  to  the  overall  cost. 

Iu  the  next  section,  the  effects  of  the  radome  shape  and  electrical 
properties  of  three  radome  materials  on  the  electrical  performance  of  the 
seeker  system  are  presented  <  :*.In  the  last  section,  a  final  choice  of  radome 
material  and  shape  is  made  and  presented . 


SECTION  IV 


ELECTRICAL  PERFORMANCE  STUDY 


1.  Introduction 

Electrical  performance  data  for  a  number  of  candidate  radome  designs 
were  determined  analytically  using  the  radome  analysis  computer  program 
developed  earlier  and  described  in  Reference  5.  The  radiating  properties 
of  the  phase  monopulse  antenna  were  characterized  using  the  plane  wave 
spectrum  representation  of  the  near  field  of  the  antenna  as  described  in 
2c  below.  The  important  electrical  performance  parameters  of  the  seeker 
system  (boresight  error,  error  slc^e,  null  depth,  and  power  loss)  were 
computed  at  three  frequencies  (8,  12,  and  18  GHz)  for  three  different 
radome  materials  (fused  silica,  polyimide,*  and  Pyroceram^  and  for  a 
number  of  fineness  ratios  (0.5  to  2,8)  and  antenna  gimbal  positions.  Graphs 
of  the  electrical  performance  parameters  of  the  antenna /radome  combinations 
are  presented  in  subsection  3  below  and  used  in  selecting  an  optimum  radome 
design  for  the  seeker  system  considered.  The  antenna/radome  geometry, 
radome  shape,  electrical  performance  parameters,  and  antenna  characteriza¬ 
tion  are  described  in  subsection  2  below. 

2.  Definitions  of  Input  and  Output  Parameters 

(a)  Radome  Shape 

The  radome  shape  considered  is  shown  in  Figure  11.  The  radome 
consists  of  a  5-inch  outside  diameter  cylindrical  section  surmounted  by  a 


The  polyimide  referred  to  in  Section  IV  is  polyimide  E-glass  which  has  a 
relative  dielectric  constant  of  3.8.  PBI  also  has  a  relative  dielectric 
constant  of  3.8;  thus  the  polyimide  curves  are  representative  of  the  PBI 
materials.  Polyimide  quartz  has  dielectric  properties  similar  to  slip- 
cast  fused  silica  (p  =  1.95  g/cc). 
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NOTE:  ANTENNA  APERTURE  DIAMETER  IS  RELATED  TO  THE  DISTANCE  R  BY  THE 
RELATION  SHOWN.  IF  RA»1.20",  THEN  1-3.70".  A 


Figure  11.  Antenna/Radotne  Geometry. 


tangent  ogive.  Various  degrees  of  streamlining  of  the  radome  are  obtained 
by  varying  the  fineness  ratio  of  the  tangent  ogive  portion  of  the  two- 
section  radome  where  the  fineness  ratio  V  is  defined  as  the  ratio  of  the 
length  L  of  the  tangent  ogive  to  its  outside  diameter  D  at  a  point  P  where 
a  tangent  to  the  ogive  shape  is  parallel  to  the  axis  of  symmetry  of  the 
radome;  i.e,,  F  ®  L/D.  In  Figure  11,  the  center  of  generation  of  the 
tangent  ogive  is  located  at  the  gimbal  point  of  the  antenna  so  that  such 
a  point  P  lies  on  the  outer  surface  of  the  radome  in  the  plane  perpendicular 
to  the  radome  axis  and  containing  the  gimbal  point.  In  the  analysis,  fine¬ 
ness  ratios  of  0.5  (a  hemisphere),  1.0,  1.5,  2.0,  2.5,  and  2.8  were  considered. 

The  location  and  orientation  of  the  near-field  aperture  of  the  circu¬ 
larly  polarized,  phase  monopulse  antenna  are  also  shown  in  Figure  11.  The 
aperture  is  offset  from  the  gimbal  point  by  the  distance  as  shown.  A 
value  of  Ra  =*  1.20  inches  was  used  in  the  analysis  to  conform  to  the  dimen¬ 
sions  of  an  actual  antenna  system  under  consideration  for  this  application. 

The  maximum  diameter  of  the  circular  aperture  which  can  fit  inside  the 
radome  is  given  by 


where 

D  «*  Outside  diameter  of  the  radome 

t  =  Thickness  of  the  radome  wall 
w 

R  =■  Distance  of  the  aperture  from  the  gimbal  point 
A 

The  maximum  wall  thicknesses  to  consider  ."ollow  from  the  equation  defining 
a  n-order  monolithic  radome  wall  [6^5  i.e.. 
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where 


n  =  1  for  half-wave  wall 
\o  =  Free-space  wavelength 

«  Dielectric  constant  of  wall  material 
9^  =  Design  angle  (60  degrees) 

Using  the  dielectric  constants  of  the  three  wall  materials  considered  in 
Equation  (2)  at  the  (center)  design  frequency  of  12  GHz  resulcs  in  the 
following  wall  thicknesses  for  first-order  walls: 


Material 

€r 

tan  6 

t  (inches) 
w 

Fused  Silica 

3.33 

0.001 

0.306 

Polyimide 

3.80 

0.006 

0,282 

Pyroceram 

5.20 

0.001 

0.233 

From  the  above  data  it  is  clear  that  the  material  with  the  thickest  wall 
is  fused  silica  (polyimide  quartz  has  dielectric  properties  similar  to 
fused  silica).  Consequently,  a  maximum  wall  thickness  of  0.31  inches  was 
used  in  Equation  (1)  above  to  determine  the  diameter  of  the  aperture.  It 
is  important  to  choose  an  aperture  of  maximum  diameter  as  explained  in 
subsection  2c  below.  (The  loss  tangents  for  the  materials  of  interest  are 
given  in  the  above  data  for  completeness;  tan  §  does  not  enter  into  the 
calculation  of  wall  thickness.) 

To  complete  the  specification  of  the  radome  shape,  it  is  necessary  to 
include  a  radome  base  as  shown  in  Figure  11  to  represent  the  bulkhead  where 
the  radome  is  attached  to  the  missile  body.  The  orientation  of  the  aperture 
inside  the  radome  at  the  maximum  gimbal  angle  considered  is  also  shew  in 


Figure  11. 


(b)  Antenna /Radome  Orientation 

The  coordinate  systems  used  to  orient  the  antenna  inside  the  radome 
are  shown  in  Figure  12.  These  coordinate  systems  represent  a  specialization 
of  the  system  defined  in  Reference  5  which  explains  why  .three  systems  of 
rectangular  coordinates  are  used  instead  of  just  two. 

The  antenna  coordinate  system  (x^,  y^,  z^)  is  located  in  the  reference 
system  at  (x  *  o,  y  *»  0,  z  =  R^)  and  oriented  as  shown.  The  principal 
planes  of  the  antenna  radiation  patterns  are  also  defined  in  Figure  ,12  as 
the  azimuth  and  elevation  planes.  For  reference,  the  spherical  coordinates 
of  the  origin  0.  of  the  antenna  system  in  the  reference  system  are 


0.  :  (r  =  1.20,  0  *»  0° ,  9  3  9(P).  The  angle 


9QP  to  bring  the  z 


and  z,  axes  into  coincidence.  The  z.  -  axis  also  coincides  with  the  true 
A  A 

electrical  boresight  direction  of  the  antenna. 

The  orientation  of  the  radome  coordinate  system  (x^  yR,  zR)  with 
respect  to  the  reference  system  may  be  specified  by  the  spherical  coordi¬ 
nates  of  the  origin  0  of  the  radome  system.  Note  that  the  z  -  axis 

R  K 

(and  z  )  passes  through  the  or:lgin  of  the  reference  system  which  also  coincides 
with  the  gimbal  point  inside  thejadorae.  For  reference,  0  is  located  at 
r^  =  1.20  (inches)  while  0^  and  9^  are  varied  to  produce  the  required 
scanning  of  the  antenna  inside  the  radome  to  various  look  directions. 

In  order  to  simplify  the  specification  of  the  antenna  look  direction, 
the  angles  0p  and  9^  are  defined  as  shown  in  Figure  12.  These  angles  are 
related  to  0^,  0r  according  to 

<Sr  -  <Sp  +  18CP  (3) 

er  -  180=  -  0L  (4) 

When  ns  =  0,  9.  =  0,  the  axis  of  the  radome  coincides  with  the  boresight 
P  L 

direction  of  the  antenna,  and  we  say  that  the  antenna  is  looking  straight 


F Igure  l 2 . 


Coordinate  Systems  Used  in  Radome  Analysis. 

Reference  System:  (x,y, z) 

Antenna  System:  '  (x.,vA>z,) 

AAA 


Radome  System: 


cv>vV 
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ahead.  When  0^  =  0  and  g^  >  0,  we  say  that  the  antenna  is  scanned  in  the 

azimuth  plane  (+  x  ,  z  plane);  when  9  <  0,  the  antenna  is  still  scanned 
It  A  Li 

in  the  azimuth  plane  but  in  the  opposite  direction  (-  x„,  z  plane). 

«  A 

Similar  saan.  of  the  antenna  in,  the  elevation  plane  is  accomplished  by 

setting  0  **  90°  and  varying  9T .  Any  other  plane  of  scan  may  be  specified 

by  setting  0  while  varying  9  to  produce  the  desired  scan  in  that  plane; 

P  L 

e.g. ,  when  0  **  45°,  scan,  in  the  so-called  crosstalk  pJ  <ne  may  be  effected. 

P 

It  is  advantageous  to  note  that  identical  effects  will  be  produced  in 
the  antenna  radiation  patterns  when  scanned  in  either  the  azimuth  or  eleva¬ 
tion  planes  because  of  symmetry  in  the  radome  and  symmetry  in  the  patterns 
of  the  antenna  considered.  For  example,  if  the  antenna  is  scanned  in  the 
azimuth  plane,  the  azimuth  pattern  will  be  most  affected,  especially  in 
bores ight  error.  The  same  effects  would  be  produced  for  scan  in  the  eleva¬ 
tion  plane,  with  the  roles  of  the  azimuth  and  elevation  principal  planes 
being  reversed.  Note  also  that  for  scan  in  any  given  plane,  the  magnitudes 
of  the  effects  will  be  the  same  for  as  for  -9^.  Consequently,  it  is 
necessary  to  consider  scan  of  the  antenna  through  only  positive  values  of 
9^  for  any  plane  of  scan  specified  by  0^  in  order  to  determine  the  magni¬ 
tudes  of  errors  in  the  boresight  direction,  error  slope,  and  so  on.  Further¬ 
more,  tests  conducted  using  the  radome  analysis  computer  program  showed  that 
maximum  errors  occurred  in  the  azimuth  pattern  for  scan  in  the  azimuth  plane; 
hence,  it  was  concluded  that  for  the  purposes  of  selecting  the  best  radome 
design,  scan  of  the  antenna  ir.  only  the  azimuth  plane  nee.  be  considered. 

This  finding  greatly  reduced  the  number  of  look  directions  required  to  pro¬ 


duce  representative  data,  resulting  in  substantial  savings  in  computation 


1 


(c)  Antenna  Characterization 


1 


P 


Ideally,  the  far-field  radiation  pattern  of  the  antenna  mounted 

in  place,  without  radome,  would  be  known  in  both  amplitude  and  phase  for 

two  orthogonal  components  at  each  frequency  and  scan  position  of  the 

antenna.  Equally  ideal  would  be  the  knowledge  of  the  near  field  distri- 

» 

bution  of  the  antenna  mounted  in  place,  without  radcme,  both  in  ampli¬ 
tude  and  phase  for  two  orthogonal  components  of  electric  field  at  each 
frequency  and  scan  position  of  the  antenna.  In  most  cases,  however,  this 
complete  description  is  unavailable  and  some  assumptions  and  approximations 
must  be  made.  In  the  case  at  hand,  the  antenna  to  be  used  has  not  bee' 
designed,  constructed  or  measured.  It  is  known,  however,  that  the  antenna 
to  be  used  should  be  circularly  polarized  and  configured  as  a  dual  plane 
monopulse  with  sum,  elevation  plane  difference  and  azimuth  plane  difference 
patterns.  The  pattern  of  the  antenna  should  be  frequency  independent  from  8- 
18  GKz  and  would  probably  use  spirals  as  its  main  radiating  elements.  And 
finally,  the  antenna  w^-ll  be  gimballed  and  must  fit  within  a  5-inch  outer 
diameter  radome  whose  maximum  wall  thickness  is  0.31  inch.  What  follows  is 
a  synthesis  of  a  complex  vector  near-field  distribution  which  meets  these 
specifications  under  the  assumption  of  no  near  field  variation  with  antenna 
scan. 

The  far-field  pattern  of  the  antenna  is  specified  in  the  antenna 
coordinate  system  of  Figure  12.  The  far-field  components  are  denoted 
EQ(e,0,r),  fy>(9,0,r)  and  Er(0,0,r)  where  Er(e,0,r)  =  0  in  the  far  field 
of  any  antenna  of  finite  size.  Also,  in  the  far  field,  the  dependence  of 
the  electric  field  components  on  r  is  of  the  form  1/r  so  that  on  a  sphere 
of  constant  radius  about  the  antenna,  this  term  is  a  constant  and  will  be 
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suppressed.  With  the  above  simplifications  the  far  field  of  the  antenna 


is  composed  of  E,*(8,0)  and 

9  0 

Let  us  now  define  a  complex  polarization  ratio  as  follows: 

E  (0,0) 

R(0,0)  =  E^(e,0) 

0 


(5) 


Civen  that  the  polarization  ratio  has  been  specified  for  all  angles 

(9,0)  let  us  now  calculate  the  relationship  between  E  (0,0)  abd  E  (9,0) 

x  y 

in  the  far  field  as  these  are  the  components  that  will  relate  to  the 

near-field  components  which  we  ultimately  want  to  determine. 

It  is  well  known  that  the  E  (0,0)  and  E  (0,0)  components  of  the  far- 

x  y 

field  may  be  written  in  terms  of  E  (0,0)  and  E  (0*0)  under  the  assumption 

8  0 

that  Er(e,0)=O  as  follows: 

E  (0,0)  =  E  (9,0)  cos  e  cos  0  -  E  (e,0)  sin  0 

x  0  0 

E  (0,0)  =  Eo(0.0)  cos  0  sin  0  +  E .(8,0)  cos  0 

y  9  9 

Incorporating  Equation  (5)  into  Equations  (6)  and  (7)  and  dividing 

Equation  (6)  by  Equation  (4)  yields: 


(6) 

(7) 


Ex(0,0)  ^  R(q,rt)  cos  9  cos  0  -  sin  0 
Ey(9,0)  R(G,0)  cos  0  sin  0  +  cos  0 


(8) 


Thus  it  is  seen  that  the  polarization  ratio  R(g,0)  also  determines  the  ratio 

between  E  (q,0)  and  E  (0,0)  for  all  angles  (g,0).  It  is  convenient  to  change 
x  y 

the  (9,0)  coordinate  used  in  Equation  (8)  to  normalize  navenumbers  coordinates 


(k  ,k  ), 
x  y 


Hie  normalized  wavenumbers  k  ,  k  and  k  are  defined  as  follows: 

x  y  z 


k  =  sin  0  cos  0 
x 

ky  =  sin  0  sin  0 

k  =  cos  0 
z 

2  2  0 
k  +  k  +  k  »  1 
x  y  z 


(9) 

(10) 

(11) 

(12) 
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r  ^  *  ■ 


From  these  Equations  the  following  useful  relationships  may  be  derived: 


8 in  9  *  /k^  + 

J  *  y 


cos  9  =  / 1  -  k2  -  k2 
J  x  y 


sin  0 


/? T? 

a/  X  y 


cos  0 


^k2  +  k2 
x  y 


Substituting  Equations  (13)  through  (f6)  into  Equation  (8)  yields: 

Ev(k  ,k  )  R(k  ,k  )  /l  -  k2  -  k2  '  k  -  k 

’  iVno  * '  ■  y  , _ * . 

y  *’  y’  R(k  a  >  A  .  k*  -  k2  k  -  k 

x  y  -7  x  y  y  X 

For  the  case  at  hand,  the  polarization  ratio,  R(k  ,k  ),  is  equal  to 

x  y 

or  all  angles  (k^,k^) ;  i.e.,  the  radiation  is  circularly 

polarized. 

The  power  pattern  of  an  antenna  may  be  specified  as 


P(k  .k  j-'fryy  VA-V  +  aX-V  MW  + 
x  y  11 

where  T|  is  the  intrinsic  impedance  of  free  space  in  ohms  and  *  denotes 
complex  conjugate.  At  each  point  of  the  far  field,  the  field  is  a  plane 
wave  for  which 

lc  .  E  =  0 
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Equation  (19)  may  be  expanded  and  E  (k  ,k  )  may  be  determined  as  a  function 

z  x  y 

of  E  (k  ,k  )  and  E  (k  ,k  )  as  follows: 
x  x  y  y  x  y 


E  (k  ,k  )  =  ~  fk  E  (k  ,k  )  +  k  E  (k  ,k  )1 
z\  x>  y/  k  x  xv  x*  y'  y  y  x’  y  J 


Substituting  Equation  (20)  into  Equation  (18),  the  far-field  power  pattern 

may  be  determined  as  a  function  of  E  (k  ,k  )  and  E  (k  ,k  )  as  follows: 

*  x  x  y'  y  x  y 

E*(k  ,k  )  E  (k  ,k  )  k2  E*(k  ,k  )  E  (k  ,k  )  k2 

P/k  k )  = ,  2L ■  *,.  y.  n  +  -2  +  y  *  y  y  *  _y  h  +  1 . 

P(V V  1J  L1  +  k2  J  T)  L  k2  J 

2  Z 


E*(k  ,k  )  E  (k  ,k  )  +  E*(k  ,k  •  E  (k  ,k  )  _  k  -  k 
x  x  y'  y  x  y'  y  x  >  x  x  y'  r  x 


[  ^  1  (21> 


Let  us  digress  to  show  three  relationships  for  arbitrary  complex 
numbers.  Let  c  =  a  +  jb  and  d  =*  e  +  jf  be  arbitrary  complex  numbers  where 
a,  b,  e,  and  f  are  r^al.  Then 

2  2 

c*c  •»  (a  -  jb)(a  +  jb)  »  a  +  b  ,  a  real  number  (22) 

c+c*  »*  (a  +  jb)  +  (a  -  jb)  =*  2a,  a  real  number  (23) 

e*d  +  d*c  =  (a  -  jb)(e  +  jf)  -f  (c  -  jf)(a  +  jb)  (24) 

«  ae  +  jaf  -  jbe  +  bf  +  ae  +  jbe  -  jfa  +  bf 
=  2  (ae  +  be),  a  real  number 

Using  Equation  (24)  we  see  that  the  last  term  of  Equation  (21)  is  real 
and  thus  the  whole  right  hand  side  of  Equation  (21)  is  real  as  it  should  be 
for  real  power  flow.  Substituting  Equation  (17)  into  Equation  (21)  yields 


p<kx>Y 


E  (k  ,k  )  E  (k  ,k  ) 
y  x  y  y  x  y 


<M*(kx,ky)M<kx’ V 


x  y 


+  (1  "  k2)  +  (M*(kx,ky)  +  M(kx,ky)(kx,ky))] 
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2  2  2 

where  k  *  1  -  k  -  k  from  Equation  (12)  has  been  used. 
z  x  y 

Defining  the  factor  of  Equation  (25)  in  brackets  to  be  B(k  ,k  ),  Equation 

x  y 

(25)  becomes  . 


P(k  ,k  ) 
x  y 


E  (k  ,k  )  E  (k  ,k  ) 
y'  x  y*  y  x*  y' 


•  [B'kx-ky>] 


where  from  Equations  (22)  through  (24),  B(k  ,k  )  is  seen  to  be  real. 

x  y 

P(k  ,k  )  may  be  written  as  a  product  of  conjugates  as: 
x  y 

P(kx’V  "  A*(kx’V  A(kx’V  (27) 

Thus  we  may  associate  the  factors  in  Equation  (26)  with  those  in  Equation 
(27)  as 


A  0vk 


E  (k  ,k  ) 


*  .  E  (k  ,k  ) 

r  v*  it' 


v>  *  w  *  B<vy]  [ 


Solving  for  E  (k  ,k  )  yields 
y  x  y 


A(k  ,k  ) 

WV  “  /V  7=^ 


E,,(k  ,k  )  may  then  be  found  from  Equation  (17)  as 
x  »c  y 

E  (k  ,k  )  «  M(k  ,k  )  E  (k  ,k  )  (30) 

x'  x  y  x  y  y  x  y  v  ' 

Unfortunately  there  is  no  unique  way  to  accomplish  the  separation  of 

the  power  pattern  as  indicated  in  Equation  (27).  One  way  to  accomplish 

this  however,  i3  to  let  A(k  ,k  )  be  tue  complex  far-field  component  of  the 

x  y 

Ey(x»y)  component  of  a  first-guess  near-field  distribution  which  has 
characteristics  similar  to  the  specified  antenna  except  that  it  has  only 
j}  Ey(x.y)  component.  Such  characteristics  might  include  spatial  limita¬ 
tions  and  spatial  variations  which  are  characteristic  of  antennas  in  the 
class  under  consideration.  Using  the  resultant  far-field  component 
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A(kx>k^),  and  E^  in  the  far  f4eld  may  be  determinedJErom  Equations  (29) 
and  (30).  The  near  field  of  the  antenna  may  then  be  calculated  from  the 
well  known  inverse  Fourier  transform  as  follows: 


sme(x-y) 


ft" Ye  (k  ,k  )1 

L  Xv  x  y  J 


-jk.r 


Eyn £(x-y> 


FT'1[Ex(kx’ky)] 


where  FT  is  the  inverse  Fourier  transform  operation  and  E^^x,  y)  and 

E  c(x,y)  are  the  x  and  y  components  of  the  near  field  distribution, 
ynr 

respect  Lvely. 

The  procedure  developed  above  allows  one  to  determine  E^(x*y)  and 

E  ,(x, y)  near  field  distributions  which  produce  a  far  field  which  is 
ynr 

circularly  polarized  (i.e.  for  the  case  R^,#)  *  j  )  everywhere  in  the 
far  field.  Also  the  near  field  distribution  has  spatial  characteristics 
which  meet  the  specifications  the  class  of  antennas  under  consideration. 

The  choice  of  the  first-guess  near-field  distribution  represents  a 
compromise  between  a  four-spiral  antenna  and  a  frequency  independent 
antenna.  The  resultant  is  an  expanding  four-spiral  phase  monopulse 
antenna.  The  four  spirals  are  located  at  the  comers  of  a  square  and 
touch  each  other.  The  center  of  each  spiral  moves  radially  outward  and  each 
spiral  expands  inversely  with  frequency  and  directly  with  wavelength.  Each 
spiral  has  a  diameter  of  approximately  1.23>„  at  all  operating  frequencies. 
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The  size  of  Che  spirals  was  limited  by  the  space  limitations  within  the 
radome  at  8  GHz,  the  low  end  of  the  frequency  band.  The  electric  field 
distribution  for  each  spiral  is  identical,  except  for  phasing,  and  is 
symmetric  about  the  axis  of  the  spiral.  The  field  on  each  spiral  had  a 
radial  cosine  distribution  with  zero  field  in  the  center  and  peak  field 
near  the  edge  of  the  spiral. 

After  readjustment  through  the  use  of  Equations  (29)  through  (32)  the 
resulting  near  field  distributions  produced  the  sura,  elevation  difference 
and  azimuth  difference  far-field  patterns  shown  in  Figures  13,  14,  and  15, 
respectively.  Tests  showed  that  the  cross  polarization  component  of  each 
of  these  fields  (the  opposite  sense  circular)  was  no  greater  than  30  db 
below  the  primary  component  of  these  fields  throughout  the  entire  far-field 
power  patterns.  As  the  near-field  distributions  had  constant  wavelength 
dimensions  at  all  operating  frequencies,  the  far-field  power  patterns  did 
not  change  over  the  specified  band. 


(d)  Output  Parameters 

The  parameters  chosen  to  characterize  the  electrical  performance 

of  the  antenna /radome  combination  were  boresight  errors  in  the  azimuth  and 

elevation  planes,  monopulse  error  slope  in  these  two  planes,  the  depths  of 

the  difference  pattern  nulls,  and  the  power  loss  in  the  3-db  contour  of  the 

sum  pattern.  For  a  given  set  of  input  parameters,  the  radome  analysis 

computet  program  was  used  to  calculate  the  sum  power  pattern  SUM  (fc  ,k ), 

x  y 

the  azimuth  difference  power  pattern  DAZ(k  ,k  ),  and  the  elevation  difference 

x  y 

power  pattern  DEL(k  ,k  ).  The  output  parameters  are  then  calculated  from 

x  y 

these  power  patterns  as  described  in  what  follows. 
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Figure  14.  Far  Field  Power  Pattern  of  Antenna  in 
Elevation  Difference  Mode 


The  boresight  error  in  the  azimuth  plane  is  found  by  first  forming 
the  ratio 


RAZ(kx) 


DAZ(kx,0) 

SUM(kx,0) 


(33) 


and  then  converting  to  decibels 

RAZDB(kx)  =  20  log  RAZ.  (34) 

The  square  root  is  taken  to  convert  the  power  patterns  to  voltage  patterns. 

Setting  ky  =  0  in  Equation  (33)  has  the  effect  of  restricting  attention  to 

the  azimuth  plane  defined  in  Figure  12.  The  discrete  value  of  kx  for  which 

RAX(k  )  is  a  minimum  is  then  found  by  a  point-by-point  co-*' ' ..n  of  the 

values  of  RAX(k  )  in  the  vicinity  of  k  =0.  (The  true  bor<-  •  \  cec.tion 
X  x 

is  at  k  =*  0,  k  =0:  the  effect  of  the  radone  is  to  caus  :  i 

x  y 

errors  in  the  electrical  boresight  direction.)  Since  the  va*. —  of  kx 

found  may  not  coincide  with  exact  null  position,  parabolic  interpolation 

(using  three  points  on  each  side  of  the  minimum  point)  is  used  to  find  the 

exact  null  position  in  wavenumber  coordinates;  viz.,  KXNULL.  A  similar 

procedure  is  used  with  DEL(k  ,k  )  to  find  the  null  position  KYNULL  in  the 

x  y 

elevation  plane  (kx  =  0).  The  boresight  errors  are  then  given  in  milli- 
radians  by 

BSEAZ  =  Sin”1 (KXNULL)  •  1000  (35) 

BSEEL  =  Sin"1 (KYNULL)  •  1000  (36) 

where  the  sense  of  the  error  is  described  by  the  following: 

BSEAZ  <  »  0  =  rr  plane  (Figure  12) 

BSEAZ  >  =»  $  =  0  plane  (Figure  12) 
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BSEEL  <  0  =>  0  *  plane  (Figure  12) 

IS  EEL  >  0  =>  0  =  y  plane  (Figure  12) 

The  monopulse  error  slope  in  azimuth,  ESLPAZ,  is  defined  by 

ESLPAZ  =  ~ -  j^RAZ(kx)  J  kx  «  KXNULL  (37) 


Since  the  derivative  changes  signs  on  either  side  of  the  null,  the  slope 
is  actually  found  on  both  sides  of  the  null  and  an  average  is  taken  of 
the  two  magnitudes.  In  the  actual  calculations,  RAZDB(kx)  is  used  so  that 
the  error  slope  in  volts  per  degree  is  given  by 


=  Antilog  (NULLDB)  r  _d_ 

20  log  e  L  dkx  \  ?  /J 


(38) 


where 

NULLDB  =  RAZDB  (KXNULL)  * 
e  *»  Eulers  constant 

4- 

and  where  k  ,  k  indicate  that  the  derivatives  are  taken  on  either  3ide  of 
x  x 

the  null.  The  depth  of  the  null  is  just  NULLDB.  A  similar  procedure  is 
carried  out  for  these  parameters  in  the  elevation  plane. 

The  final  output  parameter  is  the  power  loss  in  the  3-dB  contour  of  the 
sum  pattern.  This  loss  is  defined  by 


PI.OSS  =  -  10  log 


/  P3DBR  \ 
V  P3DB0  / 


where 


(39) 


P3DBR  =  Power  in  3-dB  contour  of  sum  pattern  with  radome  in  place 
P3DB0  =  Power  in  3-dB  contour  of  sum  pattern  without  radome 
The  3-dB  contour  is  the  half-power  contour  on  a  sphere  for  the  sum  power 
pattern.  The  power  calculated  is  that  received  by  a  circularly  polarized 
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antenna  when  the  seeker  antenna  transmits  the  sum  pattern.  The  parameter 
PLOSS  is  an  indication  of  the  reduction  in  gain  of  the  seeker  antenna  due 
to  inner  radome  reflection,  losses  in  the  radorae  wall,  distortions  of 
the  radiation  pattern,  and  depolarization  of  the  circularly  polarized 
waves  radiated  by  the  seeker  antenna  as  caused  by  the  radome. 

3.  Computed  Electrical  Performance 

The  results  of  the  computer  analysis  have  been  plotted  and  are  pre¬ 
sented  in  this  section.  The  plots  are  grouped  into  two  main  divisions: 

The  first  division  consists  of  plots  of  the  boresight  error  in  milliradian3, 
the  monopulse  error  slope  in  millivolts  per  degree,  the  power  loss  in  deci¬ 
bels  and  tv,>  depth  of  the  difference  null  in  decibels  of  the  far  field  pattern 
with  radome  in  place,  all  as  a  function  of  the  look  angle  of  the  antenna 
within  the  radome.  These  plots  are  repeated  for  various  candidate  radome 
materials  and  fineness  ratios  at  three  frequencies.  The  second  group  of 
plots  are  principle  plane  cuts  of  a  polyimide  radome-enclosed  antenna 
operating  in  the  sura  mode,  the  azimuth  difference  mode  and  the  elevation 
difference  mode.  These  plots  are  repeated  for  two  fineness  ratios  and  for 
three  frequencies. 

The  first  group  of  plots  are  presented  in  order  to  determine  which 
radome  material  is  electrically  superior  and  which  fineness  ratio  is 
optimum.  Far-field  patterns  were  calculated  for  fused  silica,  polyimide 
and  Pyroceram"'  radomes  with  fineness  ratios  of  0.5,  1.5,  and  2.5  at  8  GHz, 

12  GHz,  and  18  CHz.  All  plots  are  a  function  of  the  look  angle  of  the 
antenna  with  respect  to  the  radome  where  zero  degree  look  angle  corres¬ 
ponds  to  the  tip  of  the  radome.  Figure  16  is  «  graph  of  the  boresight 
error  produced  by  the  fused  silica  radome  for  look  angles  from  zero  degrees 
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to  35  degrees.  Also  on  this  figure  is  a  graph  of  the  monopulse  error  slope 
with  a  fused  silica  radome  in  place.  Both  of  these  graphs  contain  three 
sets  of  data,  one  for  a  radome  of  fineness  ratio  equal  to  0.5,  one  for  a 
radcme  of  fineness  ratio  equal  to  1.5  and  one  for  a  radome  of  fineness  ratio 
of  2,5,  all  operating  at  a  frequency  of  8  GHz.  As  expected,  the  boresight 

error,  for  a  fineness  ratio  of  0.5  (a  hemisphere)  is  2ero  and  the  error  slope, 

is  constant  for  the  various  look  angles.  The  boresight  error  is  seen  to  be 
greatest  for  a  fineness  ratio  of  1.5  and  returns  to  lov7er  values  for  a  fine¬ 
ness  ratio  of  2.5,  Graphs  to  be  display  .-d  later  for  a  polyimide  quartz 

radome  show  this  variation  in  boresight  error  in  more  detail  for  interme¬ 

diate  fineness  ratios.  The  error  slope  is  seen  to  be  a  weak  function  of 
fineness  ratio  and  look  angle  and,  as  shown  later  for  polyimide,  is  only 
severe  in  the  tip  region.  Figure  17  contains  the  graphs  for  power  loss 
and  null  depth  for  a  fuzed  silica  radome  at  8  GHz  for  fineness  ratios  of 
0,5.  I  5,  and  2.5.  The  power  loss  is  seen  to  monotonically  increase  with 
fineness  ratio  and  monotonically  decrease  with  look  angle.  The  graph  of 
null  depth  shows  that  the  null  tends  to  fill  in  almost  monotonically  with 
lock  angle.  Actually,  it  will  be  seen  later  that  the  null  depth  values 
peak  at  a  look  angle  of  approximately  20  degrees.  This  look  angle  corres¬ 
ponds  to  the  peak  of  one  of  the  two  difference  pattern  main  lobes  passing 
through  the  tip  region  of  the  radome.  The  0.5  fineness  ratio  null  depth  is 
the  lowest  and  remains  constant  with  look  angle,  as  expected.  Figure  18 
and  19  are  repetitions  of  Figures  16  and  17  but  at  a  frequency  of  12  GHz,  the 
design  frequency.  The  12  GHz  behavior  is  very  similar  to  the  8  GHz  behavior 
for  the  boresight  error,  the  monopulse  error  slope  and  the  null  depth,  but 
the  losses  for  all  fineness  ratios  at  12  GHz  are  very  low.  Figures  20  and  21 


are  repetitions  of  Figure  16  and  17  but  for  a  frequency  of  18  GHz.  Small 
differences  can  be  seen,  but  the  overall  performance  is  very  similar  to  the 
8  GHz  case.  Figures  22  through  27  are  similar  to  Figures  17  through  21  except 
that  the  material  is  polyimide  instead  of  fused  silica.  All  of  the  polyimide 
graphs  contain  data  for  six  fineness  ratios;  0.5,  1.0,  1.5,  2.0,  2.5,  and 
2.8.  Also,  the  look  angle  of  zero  degrees  is  included  as  well  as  additional 
look  angles  of  12.5  degrees,  17.5  degrees,  and,  on  some  graphs,  40  degrees. 
Figure  22  shows  that  the  bores ight  error  can  be  minimized  at  two  fineness 
ratios.  A  fineness  ratio  of  0.5  is  ideal  with  zero  boresight  error  for  all 
look  angles.  Second  best  is  a  fineness  ratio  of  2.5  since  this  fineness 
ratio  produces  the  least  boresight  error  of  the  non-hemispherical  fineness 
ratios.  The  error  slope  graph  shows  that  the  error  slope  increases  with 
fineness  ratio  in  the  tip  region.  This  increase  in  error  slope  is  probably 
due  to  the  attenuation  as  a  function  of  look  angle  which  would  lower  the 
gain  of  the  difference  pattern  in  the  null  region.  Figure  23  shows  that 
radome  losses  for  polyimide  can  become  quite  severe  in  the  tip  region  for 
large  fineness  raLios.  The  null  depth  graphs  show  a  pronounced  peak  as  a 
difference  pattern  main  lobe  passes  through  the  tip  of  the  radome  at  a  look  • 
angle  of  approximately  20  degrees.  The  remaining  figures  for  polyimide  at 
12  GHz  and  18  GUz  are  very  similar  to  those  for  fused  silica  but  are  given 
in  more  detail.  Figures  28  through  33  are  repetitions  of  Figure  16  througn 
21  except  that;  the  material  is  Pyroceram^  Pyroceran^-  is  seen  to  produce 
cha  acteristios  very  similar  to  fused  silica  and  polyimide. 

Two  conclusions  may  be  reached  concerning  the  three  materials  tested. 
First,  there  Is  no  dramatic  difference  between  materials  with  regard  to 
boresight  ^yror,  error  slope  or  null  depth.  Second,  the  losses  for  fused 
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silica  and  polyimide  are  approximately  the  same  whereas  the  loss  for  Pyro-, 
ceram  is  a  few  decibels  higher. 

Two  conclusions  regarding  the  best  fineness  ratio  can  be  reached.  First, 
che  fineness  ratio  of  0.5  is  clearly  the  best  electrically.  Second,  the  next- 
best  fineness  ratio  is  2.5,  for  as  can  be  seen  from  the  boresight  error  graphs, 
the  boresight  error  becomes  very  large  for  fineness  ratios  of  1.0  and  1.5  then 
decreases  to  a  minimum  at  2.5  and  then  increases  again  for  larger  ratios. 

The  second  group  of  plots  contains  principal  plane  patterns  for  a  polyi- 
mide  radome-enclosed  antenna.  The  patterns  are  for  various  look  angles  and 
for  two  fineness  ratios.  Figures  34  through  39  are  exceptions  to  the  rule; 
these  are  principal  plane  patterns  for  an  "air"  radome  (no  radome  at  all). 
Figure  34  is  the  azimuth  plane  slice  of  the  far-field  sum  and  azimuth  differ¬ 
ence  patterns.  The  sum  pattern  is  normalized  to  zero  decibels  with  a  peak 
at  zero  degrees  azimuth.  The  difference  pattern  is  normalized  to  the  sum 
pattern  with  resulting  main  lobes  of  approximately  -3  dB  at  ±20  degrees  in 
azimuth.  Figure  35  is  the  elevation  plane  slice  of  the  far-field  sum  and 
elevation  difference  patterns.  These  patterns  are  seen  to  be  identical  to 
the  previous  figure  for  the  azimuth  plane.  This  symmetry  exists  as  the 
antenna  is  symmetric  and  the  "air"  radome  is  symmetric.  Figures  36  and  37 
are  the  same  as  Figures  34  and  35  except  that  the  frequency  is  now  12  GHz. 

The  12  GHz  patterns  are  the  same  as  the  8  GHz  patterns  since  the  antenna  is 
frequency  independent  and  so  is  the  "air"  radome.  Figures  38  and  39  are 
repetitions  of  Figures  34  and  35  for  18  GHz  and  are  also  the  same  as  the 
8  GHz  patterns.  Figures  40  and  41  are  the  8  GHz  patterns,  Figures  42  and 
43  are  the  12  GHz  patterns  and  Figures  44-  and  45  are  the  18  GHz  patterns  for 
a  polyimide  radome  with  fineness  ratio  of  2.5  at  a  look  angle  of  zero  degjrees 
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in  azimuth  and  zero  degrees  in  elevation.  In  these  figures,  the  azimuth 
plane  patterns  are  seen  to  be  the  same  as  the  elevation  plane  patterns 
since  the  radorne  is  symmetric  for  this  look  angle. 

It  is  noted  that  the  losses  are  significant  at  8  and  18  GHz  and  low 
at  12  GHz.  Figures  46  through  51  are  the  8,  12,  and  198  GHz  patterns  for 
the  polyimide  radome  with  fineness  ratio  of  2.5  at  a  look  angle  of  17.5 
degrees  in  azimuth  and  zero  degrees  in  elevation.  Figures  46  and  47  are 
typical  of  these  patterns  and  show  the  asymmetry  of  the  two  principal 
planes  since  the  radome  is  not  symmetric  at  this  look  angle.  The  most 
significant  feature  of  these  figures  is  the  rise  of  the  difference  pattern 
null  in  the  azimuth  plane.  Figures  42  through  57  show  the  same  radome  with 
the  same  fineness  ratio  for  a  look  angle  of  40  degrees  in  azimuth  and  0 
degrees  in  elevation.  Figures  58  through  63  are  the  patterns  for  a  look 
angle  of  17.5  degrees  in  azimuth  and  17.5  degrees  in  elevation.  The  sequence 
of  figures  from  Figure  64  to  Figure  87  are  a  repetition  of  Figures  40  through 
63  except  that  the  fineness  ratio  is  0.5  (a  hemisphere).  The  significant 
feature  of  the  hemispherical  radome  is  that  no  distortion  of  the  patterns 
occur  for  any  look  angle.  The  only  feature  displayed  is  the  constant  losses 
at  8  and  18  GHz. 
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FIGURE  BCRESISHT  ERROR  Affl  MOKCPULSE  ERROR  SLOPE  VERSUS  LOOK 
ANGLE  m  FINENESS  RATIO  FOR  FUSED  SILICA  AT  8  GHZ, 
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PC',®  UBS  IN  3  DB  CONTOUR  OF  SUM  PATTERN  AND  NULL 
DEFIH  OF  FDOTUUSE  RATIO  VERSUS  LOOK  ANGLE  AW) 
FINENESS  RATIO  FOR  FUSED  SILICA  AT  12  GHZ. 
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FIGURE  20,  EORESIGHT  ERROR  AND  NOHHISE  ERROR  SLOPE  VERSUS  LOOK 
ANGLE  AND  FIIOiESS  RATIO  FOR  FUSED  SILICA  AT  18  GHZ. 
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FIGURE  El .  PCM  LOSS  IN  3  DB  CONTOUR  OF  Sill  PATTERN  AND  NULL 
DEPTH  OF  IMPULSE  RATIO  VERSUS  LOOK  ANGLE  AND 
FINENESS  RATIO  FOR  FUSED  SILICA  AT  18  GHZ. 
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BORESIGHT  ERROR  AND  IMPULSE  ERROR  SLOPE  VERSUS  LOOK 
ANGLE  AND  FINENESS  RATIO  FOR  POLYIMIDE  AT  I?.  GHZ. 
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FIGURE  25.  POWER  LOSS  1H  3  DB  OOMTOUR  OF  SUM  PATTERN  AID  NULL 
DEPTH  OF  mPULSE  RATIO  VERSUS  LOOK  ANGLE  AND 
FINENESS  RATIO  FOR  POLYIMIDE  AT  \Z  GHZ.  ■ 
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FIGURE  21.  POWER  LOSS  IN  3  DB  CONTOUR  OF  SUM  PATTERN  AND  NULL 
DEPTH  OF  MONOPULSE  RATIO  VERSUS  LOOK  ANGLE  AND 
FINEESS  RATIO  FOR  PYROCERAM  AT  8  GHZ. 
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FIGURE  31.  PORER  LOSS  IN  3  DB  COOT  OF  SUM  PATTERN 
DEPTH  CF  IMPULSE  RATIO  VERSOS  LOOK  AMI 
FINENESS  RATIO  FOR  PYROCERffl  AT  12  GHZ. 
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FIGURE  at.  SUM  PND  PZIMUTH  DIFFERENCE 
PATTERNS  FOR  FINENESS  RRTIO  GF  2-5  PT 
12  GHZ  FOR  LOOK  ANGLE  GF  0.0  DEG  IN 
PZIMUTH  PND  0.0  DEG  IN  ELEVRTI GN  . 


FIGURE  37.  SuM  PNO  E,EVRTIGN  DIFFERENCE 
PATTERNS  FOR  FINENESS  RRTIO  GF  2=5  PT 
12  GHZ  FOR  LOOK  PUGLE  0F  Q.O  DEG  IN 
PZIMUTH  PND  0-0  GEG  IN  ELEVATION, 
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CO  Cl. 


F  !  GliRE  4o.  SON  PND  PZINUTH  DIFFERENCE 
PATTER S  F0R  FINENESS  RATIO  OF  2.5  RT 
8  >.•  -OR  LOOK  ANGLE  OF  0.0  DEG  IN 

PZ-i"L’-'  -NO  0.0  OEG  IN  ELEVATION. 
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BJN  PND  ELEVATION  DIFFERENCE 
-OR  FINENESS  RPTIO  OF  2-5  PT 
OR  LOOK  PNGLE  OF  C-0  OEG  IN 
iND  0.0  DEG  IN  ELEVRT I  ON  . 
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FIGURE  42.  SUM  PNO  PZMUTH  DIFFERENCE 
PRTTERNS  F6R  FINENESS  RfiT  i  0  6F  2.5  PT 
! 2  GHZ  FOR  LOOK  PNGLE  0F  0.0  DEG  IN 
fiZJMUTH  RND  0.0  DEG  IN  ELEVPT I  ON . 
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FIGURE  4G.  SUM  PNO  PZIMUTH  DIFFERENCE 
PRTTERNS  FOR  FINENESS  RRTI0  0F  2-5  PT 
8  GHZ  FOR  LOOK  ANGLE  OF  17 .5  DEG  IN 
PZIMUTH  PND  0.0  DEG  IN  ELE VRT I  ON  « 
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OR  FINENESS  RPTIO  OF  2,5  PT 
LOOK  ANGLE  OF  17.5  DEG  IN 
D  0.0  DEG  IN  ELEVATION. 
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Ai  11*11 'if 

:  I  CURE  4-8.  SUM  PND  AZIMUTH  DIFFERENCE 
PATTERNS  F0R  FINENESS  RATIO  OF  2-5  AT 
! 2  GHZ  FOR  LOOK  ANGLE  OF  17.5  DEG  IN 
AZIMUTH  “MO  0.0  DEG  IN  ELEVATION.' 
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I  DURE  So-  SUM  AND  AZIMUTH  DIFFERENCE 
PATTERNS  F0R  FINENESS  RAT  1 0  0F  2-5  AT 
18  GHZ  FQR  LOOK  ANGLE  0F  P.5  DEG  IN 
AZIMUTH  AND  0-0  DEG  IN  ELEVATION- 


FIGURE  52.  SUM  PND  AZIMUTH  DIFFERENCE 
PATTERNS  FDR  FINENESS  RATIO  OF  2.5  AT 
8  GHZ  FOR  LOOK  ANGLE  OF  40.0  DEG  IN 
AZIMUTH  PND  0.0  DEG  IN  ELEVATION. 
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FIGURE  54-.  SUM  RNO  AZIMUTH  DIFFERENCE 
PATTERNS  FOR  FINENESS  RATIO  OF  2.5  AT 
12  DiiZ  FOR  LOOK  ANGLE  OF  40.0  DEC-  IN 
AZIMUTH  AND  0-0  DEG  IN  ELEVATION . 
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i  GUKi-  dco.  SUM  PND  PZiMUTH  DIFFERENCE 
E-US  FOR  FINENESS  RPTiG  GF  2.5 
GHZ  FGR  LGCK  PNGLE  GF  40  «G  DEG 
PZIMUTH  P.ND  0.0 
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FIGURE  58.  SUM  PND  AZIMUTH  DIFFERENCE 
PATTERNS  FOR  FINENESS  RATIO  OF  2.5  A 
8  GHZ'  FOR  LOOK  ANGLE  OF  17.5  DEG  IN 
AZIMUTH  AND  45.0  DEG  IN  ELEVATION. 
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FIGURE  6,0 .  3u*1  PNO  AZIMUTH  DIFFERENCE 
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‘■UR t  64.  Sum  PND  PZIMUTH  DIFFERENCE 
PATTERNS  FQR  FINENESS  RRTIG  GF  0.5  PT 
3  GHZ  F-3R  LOOK  PNGLE  GF  0-0  DEG  IN 
P  Z i MU  T H  ~ND  0.0  DEG  IN  ELEVRT I  ON « 


ELEVPTIGN  DIFFERENCE 
ENES3  RPTIG  GF  0-5  P 
PNGLE  GF  0.0  DEG  IN 
D EG  IN  ELEVP.TION. 
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i CURE  <0(0.  SUM  AND  AZIMUTH  DIFFERENCE 
PATTERNS  FOR  FINENESS  RATIO  OF  0-5  AT 
12  GHZ  FOR  ‘LOOK  ANGLE  OF  0-0  DEG  IN 
AZIMUTH  AND  C«0  DEG  IN  ELEVATION. 
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FIGURE  78.  SUM  P.ND  PZIMUTH  DIFFERENCE 
PATTERNS  FQR  FINENESS  RATI 0  QF  0-5  AT 
12  GHZ  FOR  LQQK  PNGLE  GF  40.0  DEG  IN 
PZIMUTH  AND  0-0  DEG  IN  ELEVATION. 
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FIGURE  80.  SUM  FIND  AZIMUTH  DIFFERENCE 
PATTERNS  FOR  FINENESS  RATIO  OF  0.5  AT 
IS  GHZ  FOR  LOOK  ANGLE  OF  40-0  DEG  IN 
AZIMUTH  AMD  0.0  DEG  IN  ELEVATION. 
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FIGURE  82.  SUM  RNO  RZIMUTH  DIFFERENCE 
PRTTERNS  F0R  FINENESS  RATIO  0F  0.5  RT 
S  GHZ  FOR  LOOK  ANGLE  0F  17.5  DED  IN 
AZIMUTH  RND  45.0  DEG  IN  ELEVRT 1 0N . 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  Conclusions 

Based  on  the  data  presented  In  Sections  II  through  IV,  the  following 
conclusions  may  be  drawn: 

(a)  .Radome  Materials 

All  of  the  four  radome  materials  considered  (fused  silica, 
polyimide  quartr,  polyimide  E-glass,  and  Pyroceram)  were  found  to  have 
mechanical,  thermal,  and  electrical  properties  commensurate  with  the  mission 
of  the  radome  under  consideration.  Therefore,  the  selection  of  a  single 
material  for  this  radome  application  must  be  based  oh  small  differences  in 
electrical  performance  and  other  considerations  such  as  rain  erosion  resis- 
tance,  ease  of  fabrication,  cost,  and  handling  in  the  field. 

(b)  Radome  Shape 

Of  the  six  tangent  ogive  radome  shapes  considered,  the  two  shapes 
with  fineness  ratios  of  F  »  0.5  and  F  =>  2.5  were  found  to  be  electrically 
superior.  Aerodynamically,  only  those  radome  shapes  having  fineness 
ratios  greater  than  approximately  1.5  will  fulfill  the  minimum  range  and 
velocity  requirements.  All  radome  shapes  having  fineness  ratios  in  the 
range  0.5  to  2.8  were  found  to  be  acceptable  from  a  thermal  and  mechanical 
point  of  view. 

(c)  Wall  Structure 

Due  to  the  stringent  requirements  on  mechanical  strength  and 
maximum  radome  wall  thickness,  the  monolithic  first-order  (half-wave)  wall 
structure  was  selected. 
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2.  Recommendations 

(a)  Kadome  Composition 

It  is  recommended  that  the  radome  be  a  monolithic  first-order 
(half-wave)  wall  structure  constructed  of  polyimide  quartz  material 
having  s  dielectric  constant  of  3,40,  loss  tangent  less  than  0.006,  and 
a  thickness  of  0.299  inch.  The  thickness  is  based  on  Equation  (2)  when  a 
design  frequency  of  12  GHz  and  a  design  angle  of  60  degrees  are  used.  , 

Deciding  factors  in  this  recommendation  are  superior  electrical  perfor¬ 
mance,  ease  of  fabrication,  and  ease  of  handling  in  the  field, 

(b)  Radome  Shape 

It  is  recommended  that  the  radome  shape  be  a  tangent  ogive  sur¬ 
mounting  a  5-inch  outside,  diameter  cylinder  as  shown  in  Figire  11  with  an 
inside  diameter  of  4.402  inches.  The  fineness  ratio  of  the  tangent  ogive 
section  is  recommended  to  be  2.5  so  that  both  the  electrical  and  aerodynamic 
performances  will  be  acceptable;  however,  if  the  minimum  range  were  reduced, 
a  fineness  ratio  of  0.5  would  be  recommended  because  of  its  vastly  superior 
electrical  performance. 


REFERENCES 


Van  Driest,  E.  R. ,  "The  Problem  of  Aerodynamic  Heating,"  Aeronautical 
Engineering  Review,  October  1956. 

P»yc,  E.  E. ,  "Newtonian  Aerodynamics  for  Tangent  Ogive  Bodies  of 
revolution,"  Report  X-671-65-244,  Goddard  Space  Flight  Center,  June 
1965. 


Bomar,  S.  H, ,  et.al.,  Radcme  Materials  Evaluation,  Final  Technical 
Repo:,  t  AFAL-TR-73-222  by  Engineering  Experiment  Statior,  Georgia 
Institute  of  Technology  for  Air  Force  Avionics  Laboratory  under 
contract  F336L.-72-C-1615,  July  1973. 

Chase,  V.,  Manufacturing  Technology  for  Polybezimidazole  Resin  in 
Composite  Structures .  Technical  Report  AFML-TR-70-48,  Vol.  II,  by 
'•rsittaker  Corporation  for  Air  Force  Materials  Laboratory  under 
tract  Al -33(615)1180,  January  1971. 


Joy,  E.  B.  and  Huddleston,  G.  K. ,  "Radome  Effects  on  Ground  Mapping 
Radar,"  Contract  DAAH01-72-C-0598,  Georgia  Institute  of  Technology, 


Atlanta,  Georgia,  March  1973  and  August  1973. 


Ualt-nn  7  H  Parlnmo  WsmilhnriV  an  A  nl  ao 


Marcel  and  Decker,  Inc.;  New  York,  1970. 


,y 


